A note on color neutrality in NJL-type models 
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By referring to the underlying physics behind the color charge neutrality condition in quark 
matter, we discuss how this condition should be properly imposed in NJL-type models in a phe- 
nomenologically meaningful way. In particular, we show that the standard assumption regarding the 
use of two color chemical potentials, chosen in a very special way, is not justified in general. When 
used uncritically, such an approach leads to wrong or unphysical conclusions. 



in . 
o ■ 
o ■ 

<n : 

<u : 
m . 

(N : 

<n : 
> : 

o\ ■ 

oo : 
o . 

in ■ 
o ■ 



i 

CD . 
43 ■ 



X 



It is generally expected that cold dense quark matter 
is a color superconductor. (For reviews on color super- 
conductivity, see Ref. At asymptotically high den- 
sities, when the QCD running coupling constant is very 
small, this can be shown from first principles using weak- 
coupling techniques 2\. At not-so-large densities, how- 
ever, as present, e.g., in the cores of compact stars, the 
weak-coupling expansion breaks down. In order to get at 
least some schematic picture of the phase structure in this 
regime, one therefore has to invoke models. In this con- 
text, NJL-type models play a prominent role. (For a re- 
cent review, see Ref. 0) These models are characterized 
by the fact that the gluon fields of QCD are replaced by 
point-like quark-quark or quark-antiquark interactions. 

In general, the interactions are chosen to be consis- 
tent with the most important symmetries of QCD. How- 
ever, due to the lack of gauge fields, the local SU(3)-co\ot 
(SU(3) C ) symmetry of QCD becomes replaced by a global 
one in the NJL model. This simplification has important 
consequences: 

In QCD, the homogeneous color superconducting 
ground state is automatically color neutral. According 
to Gauss' law, any sample of color-charged matter would 
create a color-electric field, leading to a diverging en- 
ergy density in the infinite volume limit. Thus, matter 
is forced to remain color-neutral by the generation of a 
color-electrostatic potential, i.e., a non-zero gluon con- 
densate (A° n ) 7^ in one or more of the eight color com- 
ponents a Q| . (Note, however, that here and in the entire 
article we assume that the matter remains homogeneous, 
while we cannot exclude the emergence of locally colored 
inhomogeneous solutions, like mixed phases [j| or LOFF 
phases |g, as long as they are globally color neutral.) 

More recently, the above arguments have been con- 
firmed by calculating gluon condensates within the 
framework of QCD at asymptotic densities 0,0,0. For 
instance, for the standard set-up of a two-flavor color su- 
perconductor (see Eq. (JJJ below), one finds (Ag) ^ 
and (A®) — in all other components. It is obvious that 



such a condensate is well defined only after gauge fixing 
and breaks the remaining global color symmetry spon- 
taneously. In particular, this implies that it could be ro- 
tated into other color directions by global gauge trans- 
formations. 1 

In NJL-type models, on the other hand, there are no 
gauge fields which can realize the color neutrality dy- 
namically. As a consequence, the color-superconducting 
ground state is not automatically color neutral. In order 
to get a physically meaningful result , on e therefore has to 
impose color neutrality by hand Q, [HJ HH] • To this end, 
one introduces a set of color chemical potentials {^ a } 
which are chosen in such a way that the corresponding 
color charge densities vanish: 



Here 



n a = {^T a f) 



an 
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is a quark field with three color components ( "red" , 
"green", "blue"), T a = 4p denote the generators of 
SU(3) C , and Q = —y\nZ is the grand canonical po- 
tential of the system per volume. 

In the light of the above discussion, the color chemi- 
cal potentials should be identified with the correspond- 
ing components of the gluon condensate in QCD, fi a = 
g (A®). In fact, from the general path-integral expression 
of the partition function one obtains 
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Tr[S(x,x)r° a ] 



(3) 



where S is the dressed (Euclidean) Nambu-Gorkov quark 
propagator and 
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Employing a local gauge transformation, it could even be elim- 
inated completely. However, in the corresponding gauge the di- 
quark condensate would pick up a space-time dependent color 
phase, making it rather inconvenient for practical purposes. 
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is the quark-gluon vertex in Nambu-Gorkov space. Thus, 
the density n a takes the form of a tadpole, corresponding 
to a gluon field A® attached to a quark loop. 

From this perspective it appears natural to consider 
eight different color chemical potentials, a — 1, . . . , 8, and 
to impose Eq. for the corresponding eight derivatives. 
This is at variance with the "standard definition" of color 
neutrality in NJL-type models as having equal densities 
of red, green and blue quarks, 

n r = rig = ni, ■ (5) 

Here n c = (ipltpc)- This condition translates into requir- 
ing Eq. only for the diagonal generators, a = 3 and 
a = 8, 

n 3 = ~2\ n r~ n g) — and ng = (n r + n g — 2n&) = , 

(6) 

while the non-diagonal generators are ignored. 2 

However, keeping in mind the underlying physics, it 
is not always a priori clear whether the restriction of 
Eq. to ri3 and n$ is justified. In fact, since both, the 
gluon condensates and the "tadpoles" Eq. J3J| are 

subject to gauge transformations, it is obvious that it 
cannot be physically correct under all circumstances. For 
instance, one can easily rotate 713 and ng away by an 
appropriate global SU(3) C transformation. Then, since 
such a transformation should leave the physics invari- 
ant, if the original state was not color neutral, the trans- 
formed state cannot be color neutral either, even if both, 
the new 713 and the new ng, vanish. In that case there 
must be non- vanishing non-diagonal density components 
and their disregard may lead to false conclusions. As we 
will discuss , th is is exactly what happened in two recent 
papers |12l Il3| on the ground state of color neutral two- 
flavor color superconductors (2SC). 

To this end, let us consider the specific example of 
a two-flavor system in the 2SC phase. In the standard 
ansatz for this phase, only the red and green quarks par- 
ticipate in a diquark condensate, whereas the blue quarks 
remain unpaired: 

A 2 ^(^ T c 75 T 2 r 5 ^) = , 

Ai ~(V> t C 75 t 2 TW>) =0 . (7) 

Here t 2 is a Pauli matrix acting is flavor space, while the 
T a are the generators of SU(3) C , as before. If we use a 
common chemical potential /i for all colors, the resulting 
ground state is not color neutral, but the densities of the 
quarks participating in the condensate, i.e., the red and 



2 Non-diagonal color densities and the corresponding chemical po- 
tentials hav e b een considered in the early Ginzburg-Landau anal- 
ysis in Ref . , but they got practically "forgotten" in most later 
works. 



green ones, are larger than the density of the blue quarks: 



n r = n g > njj =>- ng > . (8) 

It turns out that this is the only non- vanishing color den- 
sity for this state |9j (see also below). 

In this situation, the standard procedure is to intro- 
duce a color chemical potential /ig which is chosen in 
such a way that ng vanishes. As a result, the free energy 
gets somewhat larger, i.e., the color neutral ground state 
is formally less favored than the colored one at /ig = 0. Of 
course, this should be considered as an artifact of the NJL 
model: In QCD, the presence of long-range color forces 
would lead to a diverging energy density and therefore 
prevent the emergence of homogeneous colored solutions. 
Instead, as pointed out above, a static gluon field arises 
which exactly neutralizes the color charge. Therefore, in 
the NJL model, the color-neutral NJL-model solution 
which is obtained after introducing a non-vanishing /ig 
should be considered as the physical one. 

However, it was recently claimed in Ref. 01 that this 
solution is unstable with respect to the emergence of non- 
zero condensates Ai and A 2 . This was seemingly con- 
firmed by Ref. [13l | who found a symmetric state with 
Ai = A 2 = A3 to be the most favored "color neutral" 
solution. In both references, color neutrality was defined 
as in Eq. ©, i.e., restricted to the diagonal color charges. 

This result is easy to understand. Since the origin of 
the non-vanishing ng in Eq. (JSJ) is the fact that the blue 
quarks remain unpaired in the standard ansatz Eq. , it 
is plausible that a symmetric ansatz, where all colors are 
paired in a "democratic" way, should lead to equal den- 
sities of red, green and blue quarks without introducing 
a non-zero /ig. On the other hand, the symmetric ansatz 
can be obtained from Eq. (JZJ) by a rotation in color space. 
Thus, since the NJL model Lagrangian is invariant under 
global SU (3) transformations, the rotated ground state is 
degenerate with the colored ground state of the standard 
ansatz Eq. Q) at /ig = 0. Recalling that the latter is lower 
in free energy than the color neutral solution with non- 
zero /ig, we could naively conclude that we have found a 
new color neutral solution that is more favored than the 
standard one. 

Of course, this cannot be true: The rotated ground 
state is physically equivalent to the unrotated one and 
therefore it cannot be color neutral. This means, if n3 = 
ng = for the symmetric ansatz, its color content must 
be hidden in the non-diagonal components. 

In order to analyze this in more detail, we consider a 
continuous SU (3) c transformation that rotates the con- 
ventional diquark condensate (0, 0, A3) into (A' 1; A 2 , A 3 ) 
with A'j = A 2 . This is given by 

^ _> <p> = Uip , (9) 
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with 



vanishing n a in the unrotated state one finds: 



U = exp [i<p(\ 5 - A 7 )] 



/ 



cos2 (75 



\ sin(y2y?) 

V 7T~ 



sin(\/2tyS>) 
V2 



sin( V2ip) 

cos(v / 2 | ^) / 



(10) 



The diquark condensates are then transformed as 

(ip T C l5 T 2 T a iP) 

-> (^ T C 7 5T 2 T Q V') = (V- T C 75 r 2 U T T a U^) . (11) 

Indeed, 

£/ T T 2 [/ = cos(V2» T 2 - -1= sin(V2» (T 5 + T 7 ) , 
C/ T T 5 C7 = 

-lsin(^)T 2 + cos 2 (-^)T 5 - sin 2 (-|)T 7 , 
C/ T T 7 C7 = 



-Lsin(v^)T 2 - sin 2 (-^)T 5 + cos2 (^) T 7 , 



(12) 



which yields for Ai = A 2 = 0: 



A[ = A' 2 = -p sin(v / 2(^) A 3 , A3 = cos(V2^) A 3 . 
v2 

(13) 

In particular, this implies the symmetric case 



Ai = A' 2 



A3 
x/3 



for 008(^2^) 



y/3 



(14) 



By using the same color transformation, one could in- 
troduce the following set of two generators of the Cartan 
algebra: 

T3 = U^T 3 U 

= cos(V2ip) T 3 + -^= sm(V2ip) (T 4 + T 6 ) 
v2 

% ee rfTaU 

= \ ( 1 + 3cos(2v / 2»") T s + ^ sm 2 (V2tp) T x 



(15) 



-~^|sin(2\^0 (T 4 -T 6 ) 



The other generators T a can be transformed in the analo- 
gous way. From these relations one can immediately read 
off the transformed color densities n a = T a if>') . For 
the interesting case that 713 and n 8 are the only non- 



n[ = — sin 2 (V2(p) n 8 , 
n 3 = cos(V2ip) n 3 , 



1 1 13 

j= wa{y/2(p) n 3 - - J - sin(2V2<^) n 8 , 



V2 



1 1 /3" 

6 = sin(V2<p) n 3 + - y - sin(2V2<p) n 8 



8 = 7 fl + 3cos(2v / 2^) 



(16) 



and n 2 = n' 5 = n' 7 = 0. Note that this implies that 
the length \ft\ = n 2 a of the vector ft = (ni, . . . , rig) 
remains invariant under color rotations. 
For the symmetric case, Eq. (|14|) . one finds 



«3 , 



1 

1 

~7I ( n 3 + «s) 

--7= (n3-ns) 



(17) 



and n' 2 = n' 5 = n' 7 = n' 8 = 0. Hence, for the 2SC phase 
with vanishing fi a , in which n 3 — and only n 8 7^ in 
the unrotated state, Eq. (JJJ, we indeed find that n 3 and 
n 8 vanish in the symmetric state (as well as n' 2 , n' 5 , and 
n 7 ), but 



-n 4 



(18) 



This means that the symmetric ansatz, Eq. (|14|l . without 
color chemical potentials does not lead to a color neutral 
solution, in agreement with our general arguments. 

Note that the ground state of Ref. ^3 can a l so be ob- 
tained from the conventional configuration in Eq. JJJ by 
a color rotation as defined in Eqs. © and (|10|l with an 
infinitesimally small angle (p. If done properly, one then 
has to rotate the chemical potentials as well. Since the /i a 
transform exactly as the n a , it follows from Eq. i|l(i|) that 
infinitesimal small chemical potentials /i a , a — 1,4, 6, are 
needed (in addition to /i 8 ) in order to keep the rotated 
state color neutral. In turn, the neglect of these addi- 
tional chemical potentials induces color charges which 
have been overlooked in Ref. |l2( . Of course, when cor- 
rectly neutralized, the rotated state is exactly degenerate 
with the unrotated one. 

In conclusion, we have pointed out the importance to 
ensure Eq. Q for all color components a = 1, . . . , 8 in 
order to obtain a physically meaningful description of 
color-superconducting phases in NJL-type models. We 
have demonstrated that, when used uncritically, the stan- 
dard approach of restricting Eq. JTJ to a — 3 and 8 could 
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lead to wrong results. This is in line with the underlying 
physical picture, where the various color chemical poten- 
tials /i a correspond to the time components of a static 
gluon field, which transform as members of an octet un- 
der SU(3) C . 

The following remark is in order here. Out of the eight 
generators, T%, . . . , Tg, there are only two which are mu- 
tually commuting, namely T3 and Tg. This means that 
there are only two color charges which can be measured 
simultaneously. It is obviously convenient to choose 713 
and /is, but this is not mandatory. It is easy to show 
that the rotated operators, T3 and Tg (see Eq. lJT5)l). also 
commute with each other, [T^Tg] = 0. This means that 
the corresponding color charges, which are linear combi- 
nations of the various components of the n a in the origi- 
nal basis, are as "physical" as 713 and ng. This is similar 
to the more familiar case of angular momentum: While 
it is customary to measure J3, i.e., to choose the z-axis 
as quantization axis, one could equally well quantize the 
system along any other axis. 

Also note that for any combination of n a , one can al- 
ways find a color rotation with 

(m,..., 718)^(0,0,713,0,0,0,0,4) , (19) 

i.e., it is always possible to work in a basis where n^ and 
ri8 are the only non- vanishing components (at most). In a 
given basis, however, the restriction to n^ and n% without 
detailed analysis is in general not justified. 

Finally, it is natural to ask whether the restriction of 



Eq. to a = 3 and 8 has already led to wrong conclu- 
sions in other cases in the literature. In particular, one 
should ask whether starting from Eq. and tuning the 
value of /is is sufficient to achieve color neutrality in the 
2SC phase. To answer this question, we can directly em- 
ploy Eq. 0. For the (gapless) 2SC phase, the explicit 
form of the propagator S has been determined, see, e.g., 
Refs. [2j,[l4|. Its color structure contains only two types of 
terms, namely diag c (l, 1, 0) and diag c {0, 0, 1). Thus, only 
their traces with T% can be nontrivial, i.e., there can be 
no non-zero n a , except of ns- It can be shown in a similar 
way that only T3 and Ts tadpoles appear in the (gapless) 
CFL and uSC phases [13]. Moreover, it has been shown 
in Ref. @ that only A% and A% condense in the CFL+if 
phase. Hence, at least for the most common phases, the 
standard procedure turns out to be sufficient. It is im- 
portant, however, to verify this statement case by case. 
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